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Abstract 

It is hoped that the accuracy of a variety of lattice calculations will be improved by 
perturbatively eliminating effects proportional to the lattice spacing. In this paper, 
we apply this improvement program to the heavy quark effective theory currents 
which cause a heavy quark to decay to a light quark, and renormalize the resulting 
operators to order as- We find a small decrease in the amount that the operator 
needs to be renormalized, relative to the unimproved case. 
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1. Introduction 

It is hoped that the accuracy of a variety of lattice calculations will be improved by 
the elimination of effects proportional to powers of the lattice spacing, a [1] [2] 

. This improvement program can be implemented order by order in a and the 
strong coupling g 2 . Calculationally, it is relatively easy to perform the leading order 
in either of these expansions. For matrix elements of vector currents, both types of 
leading order corrections are thought to be in the 20 to 30% range [3], so it is at 
the same time important to include them, and reasonable to stop at this order. 

Much of the application of the improvement program to hadronic matrix 
elements has been to the order a improvement of the Wilson fermion action and 
operators [4] -[6]. Alternative methods for measuring matrix elements involving 
quarks that are heavy relative to the QCD scale have been proposed [7] [8]. The 
action and operators in this method can be thought of as discretizations of the 
heavy quark effective field theory action and operators [8] [9]. In this paper, we 
study the application of the improvement program to this action and the currents 
which cause a heavy quark to decay to a light quark. 

The outline of the paper is as follows. In the next section, we discuss improve- 
ment of the heavy quark action, and in section three, we discuss improvement of the 
aforementioned currents. As is generally the case, the improved lattice currents need 
to be related to their continuum counterparts. This is done perturbatively in g 2 
where the scale is set by the lattice spacing. Although quite different in motivation, 
this perturbative calculation is similar to the renormalization of other discretizations 
of the heavy quark effective theory currents performed in reference [10]. That 
calculation and the present one rely on the framework and results of reference [11], 
which are summarized in section four before the new perturbative results are 
presented. In section five we give our conclusions. 

2. Improvement of the Heavy Quark Action 

In this section, we argue that under two criteria of improvement, there is no need to 
modify the heavy quark action. The arguments are tree level, but this is adequate 
for the leading order improvement discussed in the introduction. 
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Roughly, the on-shell improvement condition formulated by Liischer and 
Weisz [2] is that spectral quantities such as the location of single-particle poles 
should not be corrected by terms proportional to powers of the lattice spacing. 
Examining the momentum space propagator for the heavy quark action in the 
continuum [9], 

(2.1) 



Po +«e 

and comparing it with its lattice counterpart [11], 
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(2.2) 



we see that although the propagators differ at order a, the pole is located at po = to 
all orders in a. We immediately conclude that the action has no need of improvement 
under the on-shell improvement condition, to all orders in a at tree level. 

A second argument yielding this conclusion is to compare the heavy quark 
propagator in position space in an external gauge field on the lattice to the same 
thing in the continuum. t In the continuum, the propagator from x to y is [12], 

-i5 3 (y-x)6(yo-x )Pexp-ig dtA (t,x). (2.3) 

The discretization of the action with a nearest neighbor, one-sided time derivative 
follows from this propagator [11]. On the lattice, the heavy quark propagator 
from m to n in a background field is [7] [8] , 

^ nm #(n - mo)U (n - 6)^U (n - 20) f • • • U (m)l (2.4) 

In this expression, 5 nm is the three-dimensional Kronecker 5-function, 9 (no — mo) 
is 1 if no > nio and otherwise, and Uo(m) is the lattice gauge link in the time 
direction at site m. 

If the correspondence between continuum and lattice gauge fields is [2] , 

U (n) = P expiga f dt A ((n + 1 — t)a, na), (2.5) 
Jo 

then the lattice propagator perfectly reproduces its continuum counterpart. If 
one modifies the heavy quark action (by choosing any other discretization of the 



t We thank Estia Eichten for discussions leading to this argument. 
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time derivative) the position space lattice propagator (2.4) is multiplied by a c- 
number factor depending only on no— mo but is otherwise unchanged (as long as 
the discretization does not include spatial links). With the correspondence (2.5), 
there is clearly no closer approximation to the continuum propagator (2.3) than the 
propagator (2.4) obtained from the action currently in use. 



3. Improved Heavy-Light Currents 

Since the lattice heavy quark action is not in need of improvement, we review the 
situation for the light quark action. Expressions for improved lattice heavy-light 
currents quickly follow from this discussion. 

The improved Wilson fermion action proposed by Sheikholeslami and Wohlert [6] 
is obtained from an improved action with next-to-nearest neighbor couplings [4] by 
a change of variables. The resulting action is the Wilson action plus an additional 
piece, 

(XT 

ASi = -ia^-r J2 Q(n)<Ti*v9 p i*v(n)q(n). (3.1) 

P^ u (n) is the sum of plaquettes defined in reference [5] and goes to F fliy (n) in 
the continuum limit, while a^ v = [7^,7^/2. The action has the advantage that 
it only contains nearest neighbor couplings. One can most easily obtain improved 
operators for use with this action by starting with local bilinears and making the 
replacement [3] [6], 

q(n) -> q(n) - T - ^ 7^ [Up(x)q(n + n)-U li (n- fi)h(n - p)] . (3.2) 

A similar replacement is made for q(n). 

Actually, a two-parameter family of transformations, all of which yield opera- 
tors improved to order a, can be obtained by using the equations of motion for the 
light and heavy quark. The effect of applying the equation of motion to the light 
quark field has been discussed in reference [5] and the effect of applying the equation 
of motion to the heavy quark field has been studied under the guise of temporally 
split operators [10]. In either case, the operator renormalization is changed by an 
amount which comes from the self-energy. We will not consider this generalization 
further. 



The most general heavy-light bilinear in the full continuum theory is, 



J(x) = b(x)Tq(x). (3.3) 

Here T is any Dirac matrix, and q is the light quark field. In the heavy quark 
effective theory, the corresponding operator (3.3) is [9], 

tf(x)(l O)Tq(x). (3.4) 

In a Dirac basis, the two-by-four matrix preceding T takes the form (1 0) and 
projects onto the upper two rows of T. Applying the above recipe to the local 
operator b'(n)(l O)Tq(n) we are led to consider, 

&t(n)(l 0)Tq(n) - ^ f (n)(l 0)r £ 7m [U^x)q(n + fi )-U^n- A) 1 q(n - ft) . 

(3.5) 

The first term in this expression is the usual local current used to determine 
hadronic matrix elements. The calculations necessary to renormalize this part of the 
current were mostly done in previous papers [11] [13], however it receives additional 
contributions coming from the new term in the action (3.1). The remainder of the 
current, which is naively order a, will also affect the renormalization of the current 
when it appears in loop diagrams with loop momenta of order a -1 . We now turn 
to the perturbative renormalization of the improved heavy-light current with the 
improved Wilson action of Sheikholeslami and Wohlert. 



4. Renormalization of the Improved Current 

The lattice renormalization of the heavy-light current (3.5) gives the ratio of the 
lattice operator to its counterpart in the continuum theory. We can divide the 
diagrams that contribute to this ratio into two parts. Those that give heavy and 
light quark wave function renormalization, and the 1PI vertex correction diagrams. 

The difference between wave function renormalization of the heavy quark on 
the lattice and in the continuum is g 2 /(12n 2 ) times a constant e = 4.53 computed 
in references [11] and [13]. t The corresponding light quark wave function renormal- 
ization constant was calculated in [5] . The results for the self energy-graphs which 

t Here we are using the reduced value of e which is appropriate if one fits corre- 
lation functions containing the propagator (2.4) to Ae~ 

-B(n -m )a [i ] [H] [13]. 
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r 


d\ 


d 2 




1.00 


5.46 


-7.22 


-9.21 


0.75 


5.76 


-7.23 


-8.62 


0.50 


6.30 


-7.00 


-7.80 


0.25 


7.37 


-5.72 


-6.73 


0.00 


8.79 


0.00 


-6.04 



Table 1. Previously Computed r-dependent Quantities [5] [11] [13]. 

we will need are given in terms of A^ , defined in reference [5] (when consulting 
their expressions for Ae i; note that we have taken F ooi = 1-31 [14]). 

We now turn to the vertex correction diagrams of figures 1 and 2. The 
techniques we use to evaluate these diagrams have been discussed in detail in a 
number of references [10] [11], and will not be reviewed here. 

At zero external momentum, 1(a) vanishes. As noted below equation (3.5), 
the operator has two pieces: one which is the local unimproved operator, and an 
additional piece coming from improvement. Thus the contribution of 1(b) can 
be split up into two parts which we will call the "unimproved" and "improved" 
contributions here and in the following paragraph. The "unimproved" contribution 
is [11][13], 

^(di + Gda). (4.1) 

We have defined the c-number G by GT = ^qT^q. The analytical expressions for 
di and d 2 can be found in reference [11]. We tabulate the constant A Sl and the 
constants d\ and d 2 in Table 1 for several values of the Wilson mass parameter r. 
Errors in Table 1 and 2 are at most 0(1) in the last decimal place. 

The "improved" contribution from diagram 1(b) can be combined with the 
contribution from the other three diagrams in figure 1. In order to compactly quote 
the analytic expressions for this result, we use the notation of references [15] and [16] 
for the following commonly occurring combinations: 

A 1= V sin 2 ^, 
2 

A 2 = A 4 + 4r 2 A 2 , 



A 4 = XL sin 2 Z M , 



A 5 
5 



E. 2 J/x • 2 r 
sin — sin I,,. 
M 2 



(4.2) 



r 




Ad 2 


Ad[ 


Ad' 2 


1.00 


-6.64 


-5.84 


3.43 


6.16 


0.75 


-5.30 


-3.75 


2.53 


4.94 


0.50 


-3.61 


-1.87 


1.55 


3.66 


0.25 


-1.52 


-0.45 


0.55 


2.15 


0.00 


0.00 


0.00 


0.00 


0.00 



Table 2. Changes to d\ and d 2 as a function of r. 

Additional combinations, A^ , A 2 , and A4 , are the same as above except the 
sums on \i run only from 1 to 3. Given these definitions the total "improved" 
contribution from figure 1 to the vertex renormalization is, 

,2 
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where, 



12tt 2 

-3r 2 /-rf 4 /4-A 
"iff 



(Adi + GAd 2 ), 



J 7T S 



Ado 



r 3 rdHA^ 



r d a 
J n 



.(3) 



(4.3) 



(4.4) 



We now turn to the contributions depicted in figure 2 which are due to the 
additional term in the action (3.1). In these figures, the insertion of this term is 
denoted with a cross. A tadpole diagram which vanishes has not been depicted. 
The result for the two diagrams depicted is, 

,2 
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12n 2 



(Ad[ + GAd' 2 ), 



(4.5) 



where, 



Ad[ = 



dH A 4 (3 - A] 



7T Z 



4AiA 2 



dH Aj 3) (l+r 2 A^ ; 
k 4Ap ) Aj> 3) 



2a(3)> 



(4.6) 



Combining the various results in this section, we find that the ratio of the 
lattice to continuum operators is, 



1 + 



9 Z 



12n 2 



1 1 

(di + Ad! + Adi) + (d 2 + Ad 2 + Ad' 2 )G + -e - - A El - 1 



(4.7) 



The dependance on (ia has been eliminated by setting fx = 1/a. 
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5. Conclusions 

We illustrate the use of the results of the previous section for the case of most 
interest, the current which determines the B meson decay constant, In that 
case, for reasonable values of the input parameters, the ratio of the continuum 
effective theory to full theory bilinears is numerically 0.98 [9] [13]. We need to 
multiply this by the ratio given in eq. (4.7). For T = 7075, the constant G (which 
appears in (4.7)) is —1. We take g 2 = 1.8 which is appropriate for effects arising 
from the scale n/a with 1/a = 2 GeV as argued in [17]. Taking values from tables 1 
and 2 with r = 1.00, equation (4.7) gives 1.23 and the product of the two ratios is 
1.20 (as compared to 1.28 in the unimproved case). To obtain the physical value of 
fs, one divides the lattice results for the improved current by this number. 

The operator we have renormalized is corrected both to order g 2 and to order a. 
As noted in the introduction both of these corrections are thought to be in the 20 
to 30% range. Analytically, the next perturbative corrections are proportional to 
g 2 a, g 4 , or a 2 times powers of g 2 In a [3], and there is numerical evidence that these 
further corrections are at the few per cent level for currents made from two Wilson 
fermion fields [3]. Thus it is hoped that the improved currents renormalized here 
will lead to a considerably more precise lattice determination of fs and other heavy 
quark matrix elements. 

Note Added. While preparing this manuscript we learned of unpublished work 
on this subject [18] [19] referenced in a paper on lattice measurements of several 
quantities using an improved action [20]. Martinelli and Rossi [18] argue that it 
is unnecessary to change the heavy quark action up to order a 2 . This is consistent 
with our conclusion in section two that the heavy quark action is not in need 
of improvement at any order in a. The result of Borrelli and Pittori [19] cited in 
reference [20] is consistent with our factor above if one takes g 2 = 1.0 rather than 1.8. 
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Figure Captions 

Fig. 1: The vertex correction diagrams resulting from improving the light quark 
operator 

Fig. 2: The vertex correction diagrams resulting from improving the action 
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